Energies and wave functions for a soft-core Coulomb potential 
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For the family of model soft-core Coulomb potentials represented by V{r) = — ^, with the 

(r<J + /39)? 

parameters Z > Q, /3 > 0, q > 1, it is shown analytically that the potentials and eigenvalues, E^i, 
are monotonic in each parameter. The potential envelope method is applied to obtain approximate 
analytic estimates in terms of the known exact spectra for pure power potentials. For the case g = 1, 
' the Asymptotic Iteration Method is used to find exact analytic results for the eigenvalues E^i and 

corresponding wave functions, expressed in terms of Z and /3. A proof is presented establishing 
the general concavity of the scaled electron density near the nucleus resulting from the truncated 
^ ' potentials for all q. Based on an analysis of extensive numerical calculations, it is conjectured that 

the crossing between the pair of states [{Vji), [u^i')], is given by the condition v' > + 1) and 
, £' > (£ + 3). The significance of these results for the interaction of an intense laser field with an atom 

j—i • is pointed out. Differences in the observed level-crossing effects between the soft-core potentials and 

the hydrogen atom confined inside an impenetrable sphere are discussed. 
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^ ! I. INTRODUCTION 

> : 

' The Schrodinger's time-independent equation = E"^, where the Hamiltonian H is given by (in atomic units 
[ m = h = e = 1) 

^ : 1 z 

00 ■ H=--A + V,{r), V,ir) = — — ^. (1) 

o\ ' . 

introduces the family of soft-core (truncated) Coulomb potentials, Vq, useful as model potentials in atomic physics. 
The bound states are obtained in terms of three potential parameters: the coupling Z > 0, the cutoff parameter 
I /3 > 0, and the power parameter g_> 1. The specific potentials corresponding to <? = 1 and q — 2 have been 
^ . analyzed earlier [l|, [3, H, 13, S S 0,11) 111- The potential Vi represents the potential due to a smeared charge and 
^ ' may be useful in describing mesonic atoms. The potential V2 is similar to the shape of the potential due to a 
finite nucleus and experienced by the muon in a muonic atom. Extensive applications of the soft-core Coulomb 
potential, V2, have been made through model calculations to describe the interaction of intense laser fields with 
atoms [ifl, El [13, [H, [3 EH m, [I3|- The parameter /3 can be related to the strength of the laser field, with the 
range /3 = 20 — 40 covering the experimental laser field strengths [l^l- Mehta and Patil [1] have presented analytical 
solutions for the s-state eigenvalues corresponding to the Vi potential. Patil has also discussed the analyticity of 
the scattering phase shifts for two particles interacting through the potentials Vq with q = I and q = 2. Singh et al 
^ have reported a large number of eigenvalues for the states Is to 4/ corresponding to Vi and V2 for a fixed value 
of Z; these values were obtained by the numerical solution of Eq.(l) for Z = 1; scaling laws, to be discussed here in 
section 2, extend their application to other values of Z. These authors noted that the energy-level ordering satisfied 
the condition Ei,i > E^^i , where I < i . In these formulas, v is the 'principal quantum number', defined generally 
(also for non-Coulombic potentials) as, v — I, where n is the number of radial nodes plus one. For each I value, the 
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calculated energies were found to be well represented by a Ritz type formula. Exact bound-state solutions of Vi have 
been considered earlier [1, H, d, 0] , wherein only a limited number of states with a specific choice of ^ = ... 3 have 
been treated. To our knowledge, no such study for q> 2 has been reported so far. Further, the interesting possibility 
of realizing the condition E^i = E^i^i at a common value of (3 when v > and i ^ i has not been considered as 
yet, for any Vq . In view of these observations based on the review of the previous work reported on Vq, we have 
carried out a general analysis of the characteristic features of the energies and wave functions of the complete family 
of soft-core Coulomb potentials defined by Vq as a function of all its parameters. Next, the potential envelope method 
[l8| is employed to express approximate estimates of E^g showin g an interestin g ge ometric property for all Vq. Our 
choice of the Asymptotic Iteration Method (AIM) [H, [Hi [11,121, [H, [H, [11,123, ^ enables us to present general 
new analytical results on exact bound-state solutions corresponding to Vi. It is then shown analytically that the 
electron density near the nucleus generated from all Vq is always concave. Finally, the first numerical results, on the 
crossing of energy levels in the energy spectrum corresponding to Vq, with q = I and q ~ 2 are reported. The paper is 
organized as follows: scaling and montonicity laws are established in sections 2 and 3; analytical spectral bounds are 
found by means of envelope methods in section 4; in sections 5 to 7 the Asymptotic Iteration Method is summarized 
and used for the case q = 1 to find exact analytical expressions for both eigenvalues and wave functions; in section 8 
the concavity of the scaled electron density is established for all g > 1; in section 9 we discuss the characteristics of 
the crossings of the energy levels for soft-core Coulomb potentials, and also some comparisons of these results with 
those for atoms confined inside an impenetrable sphere. 



II. SCALING 



The radial equation corresponding to ([T]) may be written 

2 \^ 2r^ {ri+/3i)^J 

where ip{0) = 0. We shall take Z to be a positive real parameter and express the general parametric dependence of 
the eigenvalues in the form E ~ E,^i(Z, (3, q), in which, as we have noted above, for a given £, the 'principal quantum 
number' v is defined by = n + ^, where n is the number of radial nodes plus one. If we make the change of variables 
r —> err in (21), where cr > is constant, multiply through by cr^, and compare eigenvalues, we immediately arrive at 
the general scaling law for this class of potentials, namely: 

E{Z,P,q) = ^EiaZ,P/a,q). (3) 

The two special cases a = 1/Z and a = (3 then yield, respectively, the special scaling laws 

E{Z, P, q) = Z^E{1, Zfi, q) = ^E{ZP, 1, q). (4) 

The parameter q is not involved because the denominator of the potential always scales like length, for every q > Q. 



III. MONOTONICITIES 



The eigenvalues E^g = E{Z, /3, q) are montone in each of the three potential parameters. In fact we shall now show: 

^<0, ->0, and -<0. (5) 

The Schrodinger operator H is bounded below. This may be shown by an application of the operator inequality 
[Hj nil > l/(4r^) which yields the general spectral bound E > minr>o[l/(8r^) + Explicit upper and 

lower bounds for all the eigenvalues may be expressed in this form with the aid of the 'potential envelope method' 
[H, [m, [H, [H, [m, 35, 36,]; this will be discussed in the next section. Thus the discrete spectrum of H may be 
characterized variationally, and from this it follows that monotonicities in the potential's dependence on the parameters 
induces the same monotonicities in the eigenvalues. We therefore prove (O by establishing, in turn, the corresponding 
monotonicities in the potential. First, by inspection, we see immediately that dV/dZ < 0. In what follows it is 
convenient to write V{r) = —Z/F{r), and to note that, if s is a potential parameter, then 

dV _ Z OF 



3 



Thus dV/ds, dF/ds, and dG/ds have the same sign, where G(r) = ln(F(r)) is given by 

G{r) = i In (r + 

We see that dG/d(3 — /3^^^/(r'? + > 0; hence we have dE/df3 > 0. Finally, we consider dG/dq. We first suppose 
that r < (3 and we define a; = r//3 < 1. In terms of this new variable we find 



dG 1 , 



1 /2:«ln(a;) 
'q \\ + x'i 



Since < x < 1, it follow that dG/dq < for r < (3. An exactly analogous argument with y = j3/r for r > (3 shows 
that dG/dq < when r > f3. Thus we conclude that dE/dq < 0. We note parenthetically that this result is in contrast 
to a generalized mean M{q) of r and [3 given, for example, by 



M{q) = 



9 \ 9 



Such a mean is known I33| to be montone increasing in q: the '2' in the denominator of the expression for M{q) makes 
the difference. 

The analysis of the monotonicity of the potential with respect to q may be used to find the limiting potential as 
q-^ oo. We have for r < (3 and x — r / (3 < 1, 



lim G{r) = lim 

q — >oo q — >oo 



Similarly, for r > (3 and y ^ f3/r < 1, wc find 



lim G{r) = lim 

q—^oo q—^oc 



Hf3) 



ln(r) 



ln(l 



ln(/3). 



ln(r) 



We conclude 



lim Vir) = Vooir) = 



f ifr </3; 
f, ifr>/3. 



(6) 



Thus, for given Z and f3, the family of potentials {V{r)}'^i has an ordered set of graphs, of which Voo{r) is the lowest: 
they intersect only at r = 0. For Z = P = 1, the family of potentials is illustrated in Fig. ([T]). 



IV. ENERGY BOUNDS BY THE 'POTENTIAL ENVELOPE METHOD' 



As we mentioned above, the operator inequality —A > l/(4r^) proved in Refs. 
energy bound that is expressed by a classical formula, namely 



[23 . ISOj l immediately yields a lower 



E > min 



1 
8^ 



V{r) 



The potential envelope method allows us to construct tighter and more specific energy formulas of this type. The 
method explores the idea that, if a given potential V{r) can be written as a smooth transformation V{r) = g{h{r)) 
of a potential h{r), for which the spectral problem is solved, then the Schrodinger spectrum generated by V{r) may 
be expressed in terms of the spectrum associated with the basis potential h{r). The method was introduced in 
1980 as a technique for the many-body problem and has subsequently been developed and used also for single-particle 
problems [3l|, [sl, l33l | , in cluding the laser-dressed potential [M], which we call V2(r) in the present article, and for 
relativistic problems [H, [sH]. If the transformation function g(h) is convex, the method yields lower energy bounds; 
conversely, when g{h) is concave, the results are upper bounds. If the basis h{r) of the transformation is a pure power 
h{r) = sgn(p) r^, then the resulting energy formula has the form 



E, 



■yi = mm 

r>0 



2r2 



Z 



(7) 
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FIG. 1: The family of potentials V{r) — —l/{r'' + 1)' for 1 < g < oo (in a.u ). The graphs are non- intersecting for r > 0, and 
are ordered, decreasing as g increases. 



where Pui{p) is a constant determined by h{r). For the class of soft-core Coulomb potentials Vq{r) we study here, 
we are able to obtain lower bounds if h{r) = — and upper bounds if h[r) — r^. With the Hydrogenic basis 
p = —1, g{h) is convex for all q > 1, and Py£(— 1) = ly. For the oscillator basis p = 2, we have P^e,{2) = 2v — {i -\- ^) 
and the transformation g{h) is concave for q 2; for larger g, in this case, the formula ([7]) yields an upper 
bound provided the critical r = f found is not too small, namely for q — {3,4,5,6}, respectively, we require f//3 > 
{0.958, 1.233, 1, 356, 1.417}. This condition arises because, even though g{h) in this case {h — r^, q > 2) is not concave 
everywhere, tangential potentials of the form a + br^, which touch in the concave region, can still provide upper 
bounds provided they do not cross Vq{r) = g{r'^). These energy bounds are in effect exploratory tools. The bounds 
obey the same scaling and monotonicity properties as the exact eigenvalues. It is perhaps geometrically interesting 
that the entire set of upper and lower approximate energy curves, expressing _E as a function of the coupling Z, 
are magnifications of a single curve, with magnification factor P^: if we write Vq{r) = Zf{r), then the parametric 
equations for the energy-bound curves are given by 



Z = P' 



r3fi(r) 



E^P^ 



1 + ./('•) 

2r^ r'^f'(r) 



(8) 



V. THE ASYMPTOTIC ITERATION METHOD 



The asymptotic iteration method (AIM) was original introduced [19| to investigate the solutions of differential equa- 
tions of the form 

y" = Xo{r)y' + s,{r)y, (' = (9) 

ar 

where Xo{r) and so{r) are C°°— differentiate functions. A key feature of this method is to note the invariant structure 
of the right-hand side of ^ under further differentiation. Indeed, if we differentiate ^ with respect to r, we obtain 

y'" = Xiv' + s^y (10) 

where Ai = Ag -f Sq + -^o Si = Sq -|- SqXq. If we find the second derivative of equation ([5]), we obtain 

= Aay' + 522/ (11) 
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where A2 = A'^ + si + AqAi and S2 = s\ + sqAi. Thus, for {n + 1)*'* and {n + 2)*'' derivative of n = 1, 2, . . . , we 
have 



and 



respectively, where 

An = A^„2 + Sn~i + AoA„_i and s„ = Sn-i + soA„_i. 

From (HH) and 1131) wc have 

A„y^"+^^ - A„_iy^"+^^ = 5„y where 5„ = A„s„_i - A„_is„ 



(12) 
(13) 
(14) 
(15) 



Clearly, from US]) if y, the solution of ([9]), is a polynomial of degree n, then (5„ = 0. Further, if (5„ = 0, then 5n' ~ 
for all n' > n. In an earlier paper [19| we proved the principal theorem of the Asymptotic Iteration Method (AIM), 
namely 



Theorem 1: Given Aq and sq C°°{a,b), the differential equation has the general solution 



y{r) = exp i J oeit)dt 



C2 + C1 / exp / (Ao(r) + 2a{T))dT dt 



if for some n> 



An A„_i 



a. 



(16) 



(17) 



VI. EXACT ENERGY EIGENVALUES 

Since the earlier work of Ciftci et al [l^, AIM has been adopted to investigate a wide range of different problems 
in relativistic and non-relativistic quantum mechanics. It should be noted that, in the process of applying AIM, 
especially to eigenvalue problems of Schrodinger-type, such as one has to overcome the following two problems 



A. Asymptotic Solution Problem 



This problem deals with the conversion of the eigenvalue problem (the absence of first derivative) to standard form 
^ suitable for the application of AIM. A general strategy to overcome this problem is to use ipir) = tpa{r)f{r) where 
TpaiT) is an asymptotic solution satisfying the boundary conditions of the given eigenvalue equation and /(r) is an 
unknown function to be determine using AIM. For we note that, as r approaches cxd, the asymptotic solution 
tpooir) of ^ satisfies the differential equation 

V'^(r)«-2i?Voo(r), (18) 

which yields 

V'oo(r) ~ e-^'. (19) 
Meanwhile, as r approaches 0, the asymptotic solution V'o(f) of ([2]) satisfies the differential equation 

-V'^'M + ^^^V'o(r)«0, (20) 
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which assumes the solution 

^o(r)«r'+i. (21) 

Consequently, we may write the exact solution of ^ as 

^(r) = r^+^e-^^'-fir). (22) 

Substituting into the Schrodinger equation we find that the radial function /(r) must satisfy the differential 
equation 



(23) 



where we denote 

= -2E. (24) 
AIM is then used to solve this second-order homogeneous differential equation for f{r). 

B. Termination Condition Problem 

This problem results when the eigenvalue problem (now in the standard form for an AIM application) fails to be 
exactly solvable. Indeed, if the eigenvalue problem has exact analytic solutions, the termination condition p7p . or 
equivalently, 

Sr,{r;E) = K{r-E)sn-i{r;E) ~ X„^i{r; E)sn{r; E) = (25) 



produces, at each iteration, an expression that is independent of r. For example, if /3 = 0, Ea. (|23p is exactly solvable 
and the termination condition (j25p yields 



Tl+l 



6niE)^Y[i-Z+{k + l)a), n = l,2,... (26) 



fc=0 

Thus the AIM condition (5„(i?) = leads to 

En = -l^, n = l,2,..., (27) 

as expected for the exact solutions of Coulomb potential. For f3 0, Eq. ([23| is not exactly solvable in general, 
however, for certain values of the potential parameters, we can obtain some analytic solutions. For example, in the 
case of q = 1, we may use AIM with 

„ / l + /2(/ + l)a 2Z \ 



to obtain a class of exact solutions given by 



Z , . 

« = 7-— -T, 29 
i + n + 1 



where n = 1, 2, 3, . . . is the iteration number used by AIM, along with the conditions on the potential parameter /3 
reported in Table 1. 
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TABLE I: Exact solutions of the radial Schrodinger equation ((2]) for g = 1. 



Conditions on [3 



-Zl3 + 1 + 2 = 

Z^{1 + 2)i3^ ~ 3Z{1 + 2){l + 2,)f3 + {21 + 3){l + 3f = 

-Z'^{1 + 3){l + 2)13^ + 6Z^{1 + 4){l + + 2)f3'^ - Z{lll^ + 50« + 54)(Z + 4)^/3 
+3(/ + 2)(2/ + 3)(Z + 4)3 =0 

Z\l + A)(l + i)(l + 2)l3* - 102-=^(; + 5)(/ + 4)(/ + 3)(/ + 2)/3^ 

+Z2(35Z^ + 300/2 _^ g23^ _|_ 720)(/ + ?>fp'^ - Z{bQl^ + 381/^ + 925/ + 720)(/ + 5)^/3 

+6(2/ + 5)(2/ + 3)(Z + 2)(/ + 5)'' = 

-Z''(l + 5)(l + 4)(l + 3)(l + 2)13'' + 15Z^{1 + 6)(/ + 5)(/ + 4)(/ + 3)(/ + 2)/3'» 
-^3(85/* + 1155/^ + 5678/2 _^ -^-^ggg^ _^ 9000)(/ + 6ff3^ 
+3^2(75/'' + 952/3 _j_ 435g^2 _^ 3522/ + 6000) (/ + Gfl3'^ 
-Z(274/'' + 3073/3 ^ 12411/^ + 21492/ + 13500)(/ + 6f[3 
+30(2/ + 5)(2/ + 3)(/ + 3)(/ + 2)(/ + 6)^ = 

Z'^il + 6)(/ + 5)(/ + 4)(/ + 3)(/ + 2)/3« - 21Zm + 7){l + 6)(/ + 5)(/ + 4)(/ + 3)(/ + 2)f3^ 
+Z^{im^ + 3430/'' + 26033/3 + 95354/^ + 167976/ + 113400)(/ + if (3^ 
-3^3(245/^ + 4592/-* + 33271/3 + 116224/^ + 195300/ + 126000)(/ + if (3^ 
+^2 (1624/5 + 28182/* + 188607/3 + 608332/^ + 945783/ + 567000) (/ + ifp'^ 
-9Z(196/5 + 3004/"' + 17753/3 + 50746/^ + 70301/ + 37800)(/ + 7ff3 
+90(2/ + 7) (2/ + 5) (2/ + 3)(/ + 3)(/ + 2)(/ + 7)** = 

-Z^(/ + 7)(/ + 6)(/ + 5)(/ + 4)(/ + 3)(/ + 2)13'' + 29.Z\l + 8){l + 7)(/ + 6)(/ + 5)(/ + 4)(/ + 3)(/ + 2)/3« 
-2^5(161/® + 4284/^5 + 46109/** + 256284/3 + 773558/^ + 1198080/ + 740880) (/ + 8)2/3^ 
+2Z*(980/'* + 25263/5 + 263144/* + 1414449/3 + 4127804/^ + 6184116/ + 3704400)(/ + 8)3/3* 
-^3(6769/^ + 165501/5 + 1632238/* + 8299620/3 + 22916602/^ + 32533488/ + 18522000)(/ + 8)*/33 
+2^2(6566/'' + 148023/5 + 1343681/* + 6287868/3 + 16004408/^ + 21011364/ + 11113200) (/ + 8)5/3^ 
-9Z{U52f + 28968/5 + 232875/* + 968195/3 + 2199048/^ + 2589112/ + 1234800)(/ + 8)®/3 
+630(2/ + 7) (2/ + 5) (2/ + 3)(/ + 4)(/ + 3)(/ + 2)(/ + 8)^^ = 

Z*(/ + 8)(/ + 7)(/ + 6)(/ + 5)(/ + 4)(/ + 3)(/ + 2)13'' - ?,QZ\l + 9)(/ + 8)(/ + 7)(/ + 6)(/ + 5)(/ + 4)(/ + 3)(/ + 2)/3^ 
+6^^(91/'^ + 3150/^ + 45472/5 + 354060/* + 1601869/3 + 4198770/^ + 5883828/ + 3386880) (/ + 13^ 
-18^5(252/^ + 8519/*^ + 120015/5 + 911495/* + 4021353/3 + 10279346/^ + 14054940/ + 7902720)(/ + 9)3/35 
+3Z*(7483/'^ + 243355/^ + 3294188/5 + 24020450/* + 101717325/3 + 249667695/^ + 328201704/ + 177811200)(/ + 9)*/3* 
-18^3(3738/'' + 114755/*^ + 1464128/5 + 10054742/* + 40105738/3 + 92835203/^ + 115350696/ + 59270400)(/ + Qfp^ 
+^2(118124/'' + 3338080/* + 39154679/5 + 247223313/* + 907897077/3 _^ 1939695507/^ + 2232161820/ + 1066867200)(/ + 9)^^/3^ 
-18^(6088/'' + 152716/*^ + 1592078/5 + 8960617/* + 29441156/3 + 56502567/^ + 58655178/ + 25401600)(/ + 9)^/3 
+2520(2/ + 9)(2/ + 7)(2/ + 5)(2/ + 3)(/ + 4)(/ + 3)(/ + 2)(/ + 9)* = 



VII. WAVE FUNCTIONS 



In this section we examine the properties of the wave functions associated with the exact eigenvalues obtained in the 
previous section. We note first, for a = and (3 = the nodeless wave function of ([2]) is given by 

^{r) = r'+i e—- (1 + '-). (30) 
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Further, for a = jj^, we have the following result. 
Theorem: For the radial Schrodinger equation if 

(/ + 3) [(6 + 30 + V(2 + 0(6 + 0] 



2{l + 2)Z 



(31) 



then the exact wave function 



Pjl + 3){2l + 3)Z + y^Pjl + 3)2(2/ + 3)Z{4{1 + 3)^ - 6(2/ + 5)Z) 
2/3(^ + 3)2 (2/ + 3) 



X 1 + a 



/3(/ + 3)(2/ + 3)Z + ^f3{l + 3)2(2/ + 3)Z(4(/ + 3)^ - p{2l + 5)Z) 
2/3(/ + 3)2(2/ + 3) 



is nodeless, while, if 



(/ + 3) [(6 + 3/) - V(2 + /)(6 + /)] 



2(/ + 2)Z 

then the exact wave function of ((32)) has exactly one node located at 



1 



(/ + 3)(2/ + 3) 



2(/ + 2)(/- V(^ + 2)(/ + 6)) 



3(/ + 2) - V(' + 2)(/ + 6) + Y^2(/ + 2)(/ + 4 + V(^ + 2)(rr6)) 



(32) 



(33) 



(34) 



Proof: Substitute /(r) = (1 + sr){l + tr) in ((23| and equate the coefficients of r's, we find that a = and s — a — t 
where t satisfies 



r - at 



(/ + 2)Z2 



(3 + 20(/ + 3)2 (3 + 2/)/3 



= 0. 



Further, f(r) now reads 



/(r) = 1 + ar + (at - r )r2 = 1 



From Descartes' rule of signs, /(r) has no real roots if 



(/ + 2)Z2 



(3 + 2/)(/ + 3)2 (3 + 2/)/3 



(;+2)z'' 



and has only one positive root if 



{l+2)Z^ Z_ 

(3+2/)(i+3)2 (3+2/)/3 



(3+2/)(i+3)2 (3+2lp 

< which holds true for /3 given by ([331 



> which holds for /3 given by ^1 



□ 



In Figures 2 and 3, we plot the un-normalized wave functions for a — and (3 given by (|3ip and p3p . respectively. 
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FIG. 2: Un-normalized wave functions ((32} for the Schrodinger equation, Eq.((2]), with soft-core Coulomb potential, q = 1, 
where a = and /3 given by (|3ip for different values of / (in a.u.). 



Further results can be obtained similarly. Indeed, for a — (jqr^^j it is straightforward to show that the solution of 
takes the form /(?-) ~ {1 + {a~ k — t)r){l + kr){l + tr) with exactly no root, one root and two roots, corresponding 
to the positive zeros of 

-Z^{1 + 3)(/ + 2)f3^ + 6Z^{1 + 4){l + + 2)f3^ - Z{llf + 50/ + 54)(/ + 4)2/3 -|- 3(/ + 2)(2/ + 3)(/ + 4)^ = 0. 



9 



25 /SO 75 100 125 150 175 
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FIG. 3: Un-normalized wave functions (|32|l for the Schrodinger equation, Eq.©, with soft-core Coulomb potential, q = 1, 
where a = 7^ and f3 given by (|33p for different values of I (in a.u.). 



VIII. CONCAVITY OF THE SCALED ELECTRON DENSITY NEAR THE NUCLEUS 

Using the radial wave function for the i^l state given by 



nui[r) — ^21+1 ' 
with the spherically averaged density normalized according to 

47r j R^i'^{r)r'^dr = An J Q{r)r'^dr = 1 , 

we define the scaled density as 



We note here that the Kato-Steiner cusp condition [33, l39|, |40|, |41| in terms of the density is given by 

9 7 



(35) 



(36) 



(37) 



(38) 



Now, for a spherical potential of the form — — + _B + /(r), where /(r) —> as r ^ it has been shown |43.l43| that 



^r(o) 



2? + 3 

Writing the soft-core Coulomb potential as 



r(4/ + 5) + 2(B-£;„ 



(39) 



we see that A = 0, S = -Z//3, and /(O) = 0. Thus we have from ((39)) 



'7r(0) = -: 



^;/(0) 



2Z + 3 

The sign of 77" (0) is then completely determined by the quantity E^ji -f ^. Writing 

1 



(40) 



(41) 



we find that 



E^i > (V-w, yg(r)Vw) > niin [Vg{r)] = - 



13 



(42) 
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Therefore, for all Vq{r), /^"(O) < i.e. the electron density near the nucleus remains concave. Specifically for Vi{r), 
the exact nodeless wave functions at /3 = ^ ^ ' correspond to i? = — 2{f+2Y '^hich leads to the simple result that 



^((0) a + 2)2- 

The corresponding ratios for the results obtained in Table I, given by AIM, can be similarly derived. Our numerical 
analysis on the trends in the eigenvalues for the potential V2 (r) , support the earlier result jl6| that the nodeless wave 



functions with energy E = — 2[i+2)'^ ^"-^ produced a.t (3 — y - ^2 ■ This leads to the ratio 

v'/{0) ^ 2Z^[V2{1 + 2) ~ 1] 
77,(0) " (2/ + 3)(? + 2)2 • 

Such ratios can be obtained for the other (3 values at which the exact energy values are known [l^ . It has been argued 
[45} that the electron density in terms of the ratio = — carries information on the identity of the nuclei, 
i.e., the value of Z and its location. The corresponding external potential at the equilibrium density thus provides 
a rationalization for the existence of the ground state energy-density functional [46^. An interesting situation arises 
in the case of the Vq{r) class of potentials when ?7;(0) vanishes. The information on the identity of the nuclei is then 

contained in the ratio We conclude this section by noting that Eq.(38) can be used to obtain the condition on 

^'^ jj^j corresponding to the critical potential at which E^i = 0. 



IX. ENERGY DEGENERACIES FOR /? > AND THE CHARACTERISTICS OF ENERGY-LEVEL 

CROSSINGS 

In this section we shall discuss the characteristic features associated with the crossings of the energy levels. Our 
results are derived from accurate numerical calculations carried out for Vq{r), q = 1 — A over /3 = — 100 in steps 
of 1. We have employed the generahzed pseudo-spectral (GPS) Legendre method with mapping, which is a fast 
algorithm that has been tested extensively and shown to yield the eigenvalues with an accuracy of twelve dig its after 
the decimal. A more detailed account, with several applications of GPS, can be found in [47l. l48l. [iol [sol. [slTIS^ and 
the references therein. In the present work, we have also verified the accuracy of these results, in a few selected cases, 
by using AIM. In order to discuss the characteristic features of the level crossings, it is useful to compare the soft-core 
Coulomb potentials with the case of hydrogen atom confined inside an impenetrable spherical cavity of radius R. 
The model potential for the spherically confined hydrogen atom , SCHA, |53| is given by VscHAir) — for r < 

R; and = 00 for r > R. The choice of parameters, R = 00 and f3 — m the potentials VscHAi''') and Vq{r) reduce 
them to that of the free hydrogen-like atom which is characterised by the well known accidental degeneracy of energy 
levels. At finite values of the parameters, the accidental degeneracy is lifted for the two potentials. The relative 
ordering of {i^tI), however, has a different functional dependence on the parameters. We note here that the level 
ordering for Vq{r) is similar to that of muonic atoms where the state with angular momentum (l + 1) is more strongly 
bound than the one with I which is different from the aufbau principle corresponding to the neutral atoms in the 
periodic table. For the two potentials VscHA{r) and Vq{r), the higher I states get relatively less destabilized as /3 or 
increases (that is to say, the eigenvalues decrease). These properties, in conjunction with the property of monotonicity, 
give rise to the crossing of a pair of states [(z^, ^) (v'jf)] with ly' > i/, V > I which are different for VscHAif) and 
Vq{r). In general, such crossings produce new degeneracy conditions not present in the free atom. In the SCHA, the 
simultaneous degeneracy condition 54 , 55i] is obtained such that all v > £ + 2, each (i^, ^) SCHA state is degenerate 
with {v + 1, i + 2) state, when both of them are confined at i? = {£ + !){£ + 2). The latter defines the radial node 
in the free {£ -j- 2,£) state. For example, at i? = 2a.u., all pairs of states (2s, Sc?), (3s,4(i), ... become simultaneously 
degenerate. For the soft-core Coulomb potential, Vq{r), our numerical analysis shows that (a) the crossing takes place 
at a certain (3 between a pair of states, [(j^, ^), {v' ,£')], defined by u' > {v + 1) and £' > {£ -\- 3), (b) (3 varies as v 
changes and it is not related to the location of radial nodes in the free hydrogen atom, (d) due to the muonic atom 
ordering, the higher than £ -I- 3 states which lie below have already crossed over before an [(z^, £),{v + l,£ + 3)] level 
crossing occurs. In Fig. 4 we have displayed the crossings of the 6s level with 7f,lg, and 7i levels given by the Vi{r) 
potential as (3 is varied. 

Similar level crossings are obtained in the cases of Vq{r) with q = 2 — 6. In Fig. 5 shows the level crossings of the 
(4p — 5.g), (5p — 6(7), (6p — 7g) and (7p, 8.g) pairs derived from V2(r) as a function of /3. This potential gives a good 
approximation for the laser-dressed hydrogen atom in intense laser field, and the parameter (3 can be identified with 



11 



Level crossings in Vi(r) potential : 6s-7(f,g,h,i) levels 

beta 

5 10 15 20 25 30 35 40 45 50 




FIG. 4: Crossings of energy levels in Vi(r) = — + P) for the 6s level with 7/, 7g, and 7i levels as a function of /3 (in a.u.). 
The energy level ordering over the range /3 = 20 — 40 is changed relative to that corresponding to lower /3 values. 

the field strength. We note in Fig. 5 that the level ordering gets inverted while passing through the range of the 
experimentally important region defined by /3 = 20 — 40. In addition to their relevance in describing the interaction of 
intense laser fields with atoms using soft-core Coulomb potential, such theoretical predictions may have applications 
in the design of optical devices. Finally, as noted in Eq.(A), the parameter q is not involved in the energy scaling. We 
have also studied numerically the variation of E^e as a function of q at fixed /? values. We do not observe any level 
crossing in these energy curves. 



-0.004 



-0.009 - 



-0.014 



LU -0.019 



-0.024 - 



-0.029 - 



Level crossings in V2(r) potential: 4p-5g,5p-6g,6p-7g pairs 

beta 

10 20 30 40 50 60 70 



-0.034 




FIG. 5: Crossings of energy levels in V2{r) — — l/(r^ + /J^)^''^ for pair of levels (4p — 5g), (5p — 6g), (6p — 7g) as a function of 
/3 (in a.u.). The energy level ordering over the range /3 = 20 — 40 undergoes changes as a consequence of cross overs. 
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X. CONCLUSIONS 

In this work we have obtained several novel analytic results on the behavior of energies and wave functions for a class 
of experimentally useful model soft-core Coulomb potentials represented by Vq{P,Z,q). Since the Hamiltonian H is 
bounded below, the energy spectrum can be characterized variationally. This allows us to conclude that monotonicities 
of potential in the parameters {Z, f3, q} induces corresponding monotonicities in the ei genv alues E^i. Upper and lower 
energy bounds have been derived by an application of the potential envelope method [18] wherein Vq{r) is written as 
a smooth transformation of a monomial potential. Such bounds are found to obey the same scaling and monotonicity 
properties as the exact E^^ and exhibit interesting geometric property described by a common magnification factor. 
The Asymptotic Iteration Method (AIM) is used for the case Vi to provide exact analytical bound-state solutions for 
both eigenvalues and wave functions. We have proved that the electron density near the nucleus generated from all 
Vq is always concave. This result, along with the specific condition on /? corresponding to the exact energy obtained 

by use of AIM for Vi, and a similar estimate, derived numerically, for V2, suggests that the ratio contains 
the information on the identity and location of the nuclear charge Z. Finally, the crossing of energy levels in the 
energy spectrum of Vq, have been analyzed using extensive numerical calculations of E,jg as a function of f3. It is 
conjectured that a pair of states, {ly'ji')], obey the crossing condition given by v' > {1/ + 1) and £'>{£ + 3). 

This condition is completely different from that of the hydrogen atom confined inside an impenetrable spherical cavity, 
where, at a characteristic common value of the radius of confinement, all such pairs of states become simultaneously 
degenerate when £' = {i + 2). Finally, in the case of the potential V2, which approximates the laser-dressed hydrogen 
atom potential, it is pointed out that the level crossing effects significantly alter the level ordering over the range of 
/3 = 20 — 40, that is identified with the field strengths of the experimental sources of intense laser radiation. 
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